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The mixed spin-1/2 and spin-3/2 Ising model on the union jack lattice is solved by establishing a mapping 
correspondence with the eight-vertex model. It is shown that the model under investigation becomes exactly 
soluble as a free-fermion eight-vertex model when the parameter of uniaxial single-ion anisotropy tends to 
infinity. Under this restriction, the critical points are characterized by critical exponents from the standard 
Ising universality class. In a certain subspace of interaction parameters, which corresponds to a coexistence 
surface between two ordered phases, the model becomes exactly soluble as a symmetric zero-field eight- 
vertex model. This surface is bounded by a line of bicritical points having interaction-dependent critical 
exponents that satisfy a weak universality hypothesis. 
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1 Introduction 

Investigation of phase transitions and critical phenomena belongs to the most intensively studied topics 
in the equilibrium statistical physics. A considerable progress in the understanding of order-disorder phe- 
nomena has been achieved by solving planar Ising models which represent valuable exceptions of exactly 
soluble lattice-statistical models with a non-trivial critical behaviour Although phase transitions of 
planar Ising models have already been understood in many respects there are still a lot of obscurities con- 
nected with a criticality of more complicated spin systems exhibiting reentrant transitions, non-universal 
critical behaviour, tricritical phenomenon, etc. It is worthy to mention, however, that several complicated 
Ising models can exactly be treated by transforming them to the solvable vertex models. A spin-1/2 Ising 
model on the union jack (centered square) lattice, which represents a first exactly soluble system exhibit- 
ing reentrant transitions 0, can be for instance reformulated as a free-fermion eight- vertex model 0. 
It should be also pointed out that an equivalence with the vertex models have already provided a precise 
confirmation of the reentrant phenomenon in the anisotropic spin-1/2 Ising models on the union jack lattice 
1 4 1, generalized Kagome lattice |5 1 and centered honeycomb lattice |6| as well. 

Despite the significant amount of effort, there are only few exactly soluble Ising models consisting 
of mixed spins of different magnitudes, which are usually called also as mixed-spin Ising models. A 
strong scientific interest focused on the mixed-spin systems arises partly on account of much richer critical 
behaviour they display compared with their single-spin counterparts and partly due to the fact that they 
represent the most simple models of ferrimagnets having a wide potential applicability in practice. Using 
the extended versions of decoration-iteration and star-triangle transformations, Fisher [7| and Yamada [8 1 
derived exact solutions of the mixed spin-1/2 and spin-S* (S > 1) Ising models on the honeycomb and 
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Fig. 1 Schematic representation of the the mixed spin-1/2 and spin-3/2 Ising model on the union jack (centered 
square) lattice. Solid (broken) lines depict the nearest-neighbour (next-nearest-neighbour) interactions. 

dice lattices, as well as on the decorated honeycomb, rectangular, triangular, and dice lattices. Notice that 
these mapping transformations were later on further generalized in order to account also for the single-ion 
anisotropy effect. The influence of uniaxial single-ion anisotropy have precisely been investigated on the 
mixed-spin honeycomb lattice 1 9 1 and on some decorated planar lattices 1 10 1, while an effect of the biaxial 
single-ion anisotropy has been explored just on the mixed-spin honeycomb lattice [ 1 1 1. To the best of our 
knowledge, these are the only mixed-spin planar Ising models with generally known exact solutions except 
several mixed-spin systems formulated on the Bethe (Cayley tree) lattices, which can be accurately treated 
within a discrete non-linear map 1 12 1 or an approach based on exact recursion equations 1131 . 

One of the most outstanding findings to appear in the phase transition theory is being a non-universal 
critical behaviour of some planar Ising models that is in obvious contradiction with the idea of universality 
hypothesis 1141 . The mixed spin-1/2 and spin-1 Ising model on the union jack lattice 1151 represents 
a very interesting system from this viewpoint as it exhibits a remarkable line of bicritical points which 
have continuously varying critical exponents obeying the weak universality hypothesis 1 16 1. Following the 
approach developed by Lipowski and Horiguchi [ 15], i.e. establishing a mapping correspondence with the 
eight-vertex model, we shall investigate in the present article the mixed spin-1/2 and spin-3/2 Ising model 
on the union jack lattice. In certain subspaces of interaction parameters, the model under investigation 
becomes exactly soluble either as a free-fermion model or a symmetric zero-field eight-vertex model. In 
the rest of the parameter space one still obtains rather reliable estimate of the criticality within so-called 
free-fermion approximation 1 17 1 when a non-validity of the free-fermion condition is simply ignored. 

The outline of this paper is as follows. In Section |2 a detailed formulation of the model system is 
presented and subsequently, the mapping correspondence that ensures an equivalence with the eight-vertex 
model will be derived. The most interesting numerical results for a critical behaviour will be presented and 
particularly discussed in Section[3] Finally, some concluding remarks are drawn in Section[4] 

2 Model system and its solution 

Let us begin by considering the mixed spin-1/2 and spin-3/2 Ising model on the union jack (centered 
square) lattice C as schematically illustrated in Fig. ^ The mixed-spin union jack lattice consists of two 
interpenetrating sub-lattices A and B that are formed by the spin-1/2 (empty circles) and spin-3/2 (filled 
circles) atoms, respectively. The total Hamiltonian defined upon the underlying lattice C reads: 

AN 2N N 

Ti-mix = — J ^ SiO-j - J' ^2 ° ko 'i - &i ' (1) 

(iJ)Cj (k.QcK i=l 
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where <jj = ±1/2 and Si = ±1/2, ±3/2 are Ising spin variables placed on the eight- and four-coordinated 
sites, J denotes the exchange interaction between nearest-neighbouring A — B spin pairs and J' labels the 
interaction between the A — A spin pairs that are next-nearest-neighbours on the union jack lattice C. 
Finally, the parameter D measures a strength of the uniaxial single-ion anisotropy acting on the spin-3/2 
sites and N denotes the total number of the spin- 1/2 sites. 

In order to obtain the exact solution, the central spin-3/2 atoms should be firstly decimated from the 
faces of sub-lattice A. After the decimation, i.e. after performing a partial trace over spin degrees of 
freedom of the spin-3/2 sites (filled circles), the partition function of the mixed-spin union jack lattice C 
can be rewritten as: 

Zmix = \\ Ufa, 17 j, (Tfc, <Jl), (2) 

{<?} i,i,k,l 

where the summation is performed over all possible spin configurations available on the sub-lattice A 
and the product is over all N faces of the sub-lattice A, which are constituted by plaquettes composed 
of a central spin-3/2 site surrounded by four spin-1/2 variables ai, aj, a^, <7; as arranged in Fig. Q The 
Boltzmann factor w(a, b, c, d) assigned to those faces can be defined as: 

w(a,b,c,d) = 2exp\j3J'(ab + bc + cd + da)/2 + /3D/4] 

(exp(2/3L>) cosh[3/3 J(a + b + c + d)/2] + cosh[/3 J(a + b + c + d)/2]}, (3) 

where (3 = l/{k&T), fee is Boltzmann's constant and T stands for the absolute temperature. At this stage, 
the model under investigation can be rather straightforwardly mapped onto the eight-vertex model on a 
dual square lattice Ct>, since the Boltzmann factor cu(a,b,c,d) is being invariant under the reversal of 
all four spin variables. Actually, there are at the utmost eight distinct spin arrangements having different 
energies (Boltzmann weights) and these can readily be related to the Boltzmann weights of the eight-vertex 
model on the dual square lattice. If, and only if, the adjacent spins are aligned opposite to each other, 
then solid lines are drawn on the edges of the dual lattice Ct>, otherwise they are drawn as broken lines. 
Diagrammatic representation of eight possible spin arrangements and their corresponding line coverings is 
shown in Fig.|2] It can be easily understood that there are eight possible line coverings around each vertex 
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Fig. 2 Eight possible spin configurations around each central spin-3/2 site and their corresponding line coverings at 
the vertices of dual square lattice. 

of the dual lattice each of them corresponding to two spin configurations, one is being obtained from the 
other by reversing all side spins. Since there is even number of solid (broken) lines incident to each vertex 
of the dual lattice £-p, the model under consideration becomes equivalent with the eight- vertex model. 
With regard to this equivalence, the partition function of the mixed-spin Ising model on the union jack 
lattice can be expressed in terms of the partition function of the eight-vertex model on the square lattice: 

Z mlx (T, J, J', D) = 2i%-»(wi,a>2, —,ug). (4) 

The factor 2 in previous equation comes from the two-to-one mapping between spin and vertex configura- 
tions (two different spin configurations correspond to one vertex configuration). 
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The Boltzmann weights, which correspond to eight possible line coverings of the eight-vertex model 
shown in Fig. [2] can readily be obtained with the aid of equation 0: 

oji = 2exp(/3D/A + /3J'/2)[exp{2/3D)cosh(3f3J) + cosh(/3J)], 

uj 2 = 2exp(/3D/4- /3J'/2)[exp(2/3D) + 1], 

u 3 = oj 4 = 2exp(/?/J/4)[exp(2/3/J) + 1], 

uj 5 = lu 6 = uj 7 = uj 8 = 2exp(/3D/4)[exp(2/3 J D) cosh(3/3J/2) + cosh(/3J/2)]. (5) 

Unfortunately, there does not exist general exact solution for the eight-vertex model with arbitrary Boltz- 
mann weights. However, if the weights (0 satisfy so-called free-fermion condition: 

uiuj 2 + = ^5^6 + ^7^8, (6) 

the eight- vertex model becomes exactly soluble as a free-fermion model treated several years ago by Fan 
and Wu 1171 . It can be readily proved that the free-fermion condition (|6} holds in our case just as D — > ±00, 
or T — > 00. The restriction to infinitely strong single-ion anisotropy consequently leads to the familiar 
phase transitions from the standard Ising universality class because of the effective reduction of the model 
system to a simple spin- 1/2 Ising model on the union jack lattice solved many years ago |2||3]|4). Within 
the manifold given by the constraint 0, the free-fermion model becomes critical as long as: 

oji + uj 2 + u} 3 + iv 4 = 2max{u/i , lu 2 , 073,0^4}. (7) 

It is noteworthy, however, that the critical condition @ yields rather reliable estimate of the criticality 
within so-called free-fermion approximation 1 17 1 even if a non-validity of the free-fermion condition l|6} is 
simply ignored. 

The second branch of exact solution occurs just as the Boltzmann weights Q satisfy the condition of 
the so-called symmetric zero-field eight-vertex (Baxter) model 1 1 1: 

wi = lj 2 , w 3 = w 4 , lu 5 =u 6 , lu 7 — lu s . (8) 

Since we already have lo 3 — U14, W5 = wg, and lu 7 = wg, hence, the symmetric case is obtained by 
imposing the condition u>i = oj 2 only, or equivalently: 

cosh(3p J) - exp(— p J') 

According to Baxter's exact solution 1 1 1, the symmetric eight-vertex model becomes critical on the mani- 
fold © if: 

oji + uj 3 + u 5 + iv 7 = 2max{wi, W3, uj 5 , uj 7 }. (10) 

It is easy to check that lo\ always represents in our case the greatest Boltzmann weight, thus, the condition 
determining the criticality can also be written in this equivalent form: 

exp(/3 c j'/2) [ exp(2/3 c L>) cosh(3/3 c J) + cosh(/? c J)] = 

1 + exp(2/3 c D) + 2exp(2/3 c /J) cosri(3/? c J/2) + 2 cosh(/3 c J/2), (11) 

where (3 C = l/(fceT c ) and T c denotes the critical temperature. It should be stressed, nevertheless, that the 
critical exponents (with exception of 5 and rj) describing a phase transition of the symmetric eight-vertex 
model depend on the function p, = 2 arctan(w 5 o;7 /ljiuj^) 1 ^ 2 , in fact: 

a = a' = 2--, /3=^-, v = u' = ^-, 7=^, 8 = 15, n = \, (12) 
/i lo/i 2/i Hp 4 

For illustrative purposes, let us explicitly evaluate the critical exponent (3 that determines disappearance of 
the spontaneous order when the critical temperature is approached from below: 

fl -i = 32 f ex P (2[3 c D) cosh(3/3 e J/2) + cosh(/3 e J/2) 1 

/J ^ arc an j [exp(2/3cZ?) + 1 ]3/4[ exp ( 2/ 3 c D) C osh(3/3 c J) + cosh(/3 c J)] 1 / 4 J ' 1 J 
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3 Results and discussion 

At first, let us turn our attention to a discussion of the most interesting results obtained for the ground- 
state and finite-temperature phase diagrams. Solid lines displayed in Fig. [3]represent ground-state phase 
boundaries separating four different long-range ordered phases that emerge in the ground state when J > 0. 




Fig. 3 Ground-state phase diagram in the J' — D plane when J > 0. Dotted rectangles schematically illustrate a 
typical spin configuration within the basic unit cell of each phase. Broken line connecting both triple points shows a 
projection of the critical line Jl H into the J' — D plane. 



Spin order drawn in dotted rectangles shows a typical spin configuration within the basic unit cell of each 
phase. As could be expected, a sufficiently strong antiferromagnetic next-nearest-neighbour interaction J' 
alters the structure of the ground state due to a competing effect with the nearest-neighbour interaction J. 
Owing to a competition between the interactions, the central spins are free to flip within the phases III and 
IV and thus, these phases exhibit a remarkable coexistence of the spin order (sub-lattice A) and disorder 
(sub-lattice E). Last but not at least, it is worthwhile to mention that a broken line connecting both triple 
points depicts a projection of the exact critical line ( II 11 into the J' — D plane. As this projection crosses 
T = plane along the ground-state transition line D/J = —3/2 — J'/2 J between the phases I and IV, it 
is quite reasonable to suspect that this line determines a location of phase transitions between these phases. 

Let us investigate more deeply this line of critical points. The critical temperatures calculated from the 
symmetric zero-field eight-vertex model must simultaneously obey both the zero-field condition (|9jl as well 
as the critical condition (^2- L 1S easv to check that the former condition necessitates —3 < J' /J < — 1 
and — 1 < D/J < 0. Fig. ^displays a projection of this critical line into the J' — T c plane (the dependence 
scaled to the left axis) and respectively, a projection into the J' — D plane which is scaled to the right axis. 
Along this critical line, the critical exponents are expected to vary with the interaction parameters as they 
have to follow the equations d!2i . For illustration, Figs. |5] [6] and [7] show how the critical indices a, (3, 
and 7, respectively, change along the critical line. Apparently, the exponents and 7 approach its smallest 
possible value 1/16 and 7/8 by reaching both triple points with zero critical temperature, while the critical 
exponent a approaches there its greatest possible value 1. It is also quite interesting to ascertain that the 
greatest values for the critical exponents (3 and 7 are slightly below the values 1/8 and 7/4, which predicts 
the universality hypothesis for planar Ising systems, while the smallest possible value of the critical index 
a is slightly above its universal value a ~ (logarithmic singularity). 

Before concluding, few remarks should be addressed to a global finite-temperature phase diagram 
plotted in Fig. [8] which displays the critical temperature as a function of the ratio J'/ J for several val- 
ues of the single-ion anisotropy D/J. Critical boundaries depicted as solid lines represent exact critical 
points obtained from the free-fermion solution under the constraint l|6} fulfilled in the limiting cases 
D/J^ ±00. Second branch of exact solution, which is related to the critical points of the symmetric 
eight-vertex model on the variety l|9}, is displayed as a rounded broken line. Dotted critical lines 
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Fig. 4 Dependence scaled to a left axis shows how 
the critical temperature changes with the ratio J'/ J, the 
curve scaled with respect to a right axis depicts variation 
of D/J along this line. Dotted lines are guides for eyes. 




Fig. 5 Dependence scaled to a left axis shows how 
the critical temperature changes with the ratio J'/ J, the 
curve scaled with respect to a right axis depicts variations 
of the critical exponent a along this line. 
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Fig. 6 The same as for Fig. [5] but the critical index (3 is Fig. 7 The same as for Fig. [5] but the critical index 7 is 
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show estimated critical temperatures calculated from the free-fermion approximation simply ignoring a 
non-validity of the free-fermion condition (|6jl for any finite value of D/ J. 

It is quite obvious from the ground-state phase diagram (Fig. [5} that a right (left) wing of the dis- 
played critical boundaries corresponds to the phase I (III) if D/J > 0, whereas it corresponds to the 
phase II (IV) if D/J < —1. Actually, the exact as well as approximate critical points resulting from the 
free-fermion solution correctly reproduce the ground-state boundaries between those phases. When the 
single-ion anisotropy term is selected within the range — 1 < D/J < (see for instance the curve for 
D/J = —0.5), however, both wings are expected to meet at a bicritical (circled) point with non-universal 
(continuously varying) critical indices as already reasoned by Lipowski and Horiguchi 1 15 1. In such a case, 
the right and left wings of critical lines separate the phases I and IV, respectively, and a line of first-order 
phase transitions is expected to terminate at this special multicritical point. There are strong indications 
supporting this concept 1 15 1, actually, the almost straight broken line depicting the zero-field condition (|9j 
should show a coexistence of these two phases as it starts from a point that determines their coexistence in 
the ground state. In addition, it is quite illuminating to see that both wings of critical temperatures referred 
to the free-fermion approximation start from this ground-state value. With regard to the aforementioned 
arguments one may conclude that a coexistence surface between the phases I and IV lies inside the area 
bounded by the line of bicritical points (rounded broken line) having the non-universal critical exponents. 
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Fig. 8 Critical temperature plotted against the ratio J'/ J for several values of D/ J. For details see the text. 

Finally, we should remark a feasible appearance of the reentrant transitions that can be observed in the 
critical boundaries nearby the coexistence points J'/ J = —3 and —1. It is quite apparent that the observed 
reentrance can be explained in terms of the coexistence of a partial order (sub-lattice A) and partial disorder 
(sub-lattice B) that emerges in both the high-temperature reentrant phases III and IV. As a matter of fact, 
the partial disorder on the sub-lattice B can compensate a loss of entropy that occurs in these phases as a 
result of a thermally induced partial ordering on the sub-lattice A in agreement with a necessary condition 
conjectured for the appearance of reentrant phase transitions pl|6). 

4 Concluding Remarks 

The work reported in the present article provides a relatively precise information on the critical behaviour 
of the mixed spin- 1/2 and spin-3/2 Ising model on the union jack lattice by establishing a mapping cor- 
respondence with the eight-vertex model. The main focus of the present work has been aimed at exam- 
ination of the criticality depending basically on the single-ion anisotropy as well as the competing next- 
nearest-neighbour interaction. The location of the critical boundaries has accurately been determined from 
corresponding solutions of the free-fermion model and the symmetric zero-field eight-vertex model, re- 
spectively, whereas the mapping correspondence is being restricted to the certain subspaces of interaction 
parameters where it holds precisely. In the rest of parameter space, the free-fermion approximation has 
been used to estimate the critical boundaries as this method should provide rather meaningful approxima- 
tion giving reliable estimate to the true transition temperatures. 

The greatest theoretical interest in the model under investigation arises due to the remarkable critical line 
consisting of bicritical points, which bounds a coexistence surface between two long-range ordered phases. 
The bicritical points can be characterized by the non-universal interaction-dependent critical exponents that 
satisfy the weak universality hypothesis. Moreover, the same arguments as those suggested by Lipowski 
and Horiguchi 1 15 1 have enabled us to identify the zero-field condition (|9jl with a location of the first-order 
transition lines separating the two ordered phases. 

It should be remarked that the considered system also shows reentrant phase transitions on account of 
the competition between the nearest- and next-nearest-neighbour interactions. Our results are in agreement 
with the conjecture 1 5 1 stating that the reentrance appears as a consequence of the coexistence of a partial 
order and disorder, namely, the partial disorder that appears on the sub-lattice B can compensate the loss 
of entropy which occurs on behalf of the partial ordering on the sub-lattice A in both the high-temperature 
partially ordered (disordered) phases. 
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